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Q^, Abstract 

1^ ' We introduce a method which allows one to recover the equations of motion of a class of 

, nonholonomic systems by finding instead an unconstrained Hamiltonian system on the full 

phase space, and to restrict the resulting canonical equations to an appropriate submanifold 
of phase space. We focus first on the Lagrangian picture of the method and deduce the 
corresponding Hamiltonian from the Legendre transformation. We illustrate the method 
with several examples and we discuss its relationship to the Pontryagin maximum principle. 
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00 ! 1 Introduction 
O 

. ^ i The direct motivation of this paper lies with some interesting results that appeared in the paper 

^ ' wherein the authors propose a way to quantize some of the well-known classical examples 

■ of nonholonomic systems. On the way to quantization, the authors propose an alternative 

Hamiltonian representation of nonholonomic mechanics. In short, the authors start off from 
the actual solutions of the nonholonomic system, and apply a sort of Hamilton-Jacobi theory 
to arrive at a Hamiltonian whose Hamilton's equations, when restricted to a certain subset of 
phase space, reproduce the nonholonomic dynamics. Needless to say, even without an explicit 
expression for the solutions one can still derive a lot of the interesting geometric features and of 
the qualitative behaviour of a nonholonomic system. However, the "Hamiltonization" method 
introduced in [17j is not generalized to systems for which the explicit solution is not readily 
available, and hence cannot be applied to those systems. 

In this paper, we wish to describe a method to Hamiltonize a class of nonholonomic systems that 
does not depend on the knowledge of the solutions of the system. Instead, we will start from 
the Lagrangian equations of motion of the system and treat the search for a Hamiltonian which 
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Hamiltonizes the dynamics as the search for a regular Lagrangian. That is, we will explain how 
one can associate to the nonholonomic equations of motion a family of systems of second-order 
ordinary differential equations and we will apply the inverse problem of the calculus of variations 
[12\ [27] on those associated systems. If an unconstrained (or free) regular Lagrangian exists for 
one of the associated systems, we can always find an associated Hamiltonian by means of the 
Legendre transformation. Since our method only makes use of the equations of motion of the 
system, it depends only on the Lagrangian and constraints of the nonholonomic system, but not 
on the knowledge of the exact solutions of the system. 

A system for which no exact solutions are known can only be integrated by means of numerical 
methods. In addition to the above mentioned application to quantization, our Hamiltonization 
method may also be useful from this point of view. A geometric integrator of a Lagrangian 
system uses a discrete Lagrangian that resembles as close as possible the continuous Lagrangian 
(see e.g. [23]). On the other hand, the succes of a nonholomic integrator (see e.g. [9l[l5]) relies 
not only on the choice of a discrete Lagrangian but also on the choice of a discrete version of the 
constraint manifold. It seems therefore reasonable that if a free Lagrangian for the nonholonomic 
system exists, the Lagrangian integrator may perform better than a nonholonomic integrator 
with badly chosen discrete constraints. Work along these lines is in progress. 

It should be remarked from the outset that the Hamiltonization we have in mind is different 
from the "Hamiltonization" used in e.g. the papers O [T3]. Roughly speaking, these authors 
first project a given nonholonomic system with symmetry to a system on a reduced space and 
then use a sort of time reparametrization to rewrite the reduced system in a Hamiltonian form 
in the new time (this is the so-called Chaplygin's reducibility trick). In contrast, we embed the 
(unreduced) nonholonomic system in a larger Hamiltonian one. 

In the second part of the paper, we show that in the cases where a regular Lagrangian (and thus 
a Hamiltonian) exists, we can also associate a first order controlled system to the nonholonomic 
system. As an interesting byproduct of the method it turns out that if one considers the optimal 
control problem of minimizing the controls for an appropriate cost function under the constraint 
of that associated first order controlled system, Pontryagin's maximum principle leads in a 
straightforward way to the associated Hamiltonians. 

We begin with a quick review of nonholonomic mechanics in section 2, where we introduce some 
of the well-known classical nonholonomic systems which fall into the class of systems we will be 
studying in the current paper. We then begin our investigations in section 3 with the Lagrangian 
approach to the problem. We detail the various ways to associate a second-order system to a 
nonholonomic system, which then forms the backbone of our subsequent analysis. In section 4 
we briefly review the set up for the inverse problem of the calculus of variations, and then apply 
it to some of the associated second-order systems. We derive the corresponding Hamiltonians 
in section 5 and discuss their relation with Pontryagin's maximum principle in section 6. At the 
end of the paper we provide a few directions for future work on generalizing our findings to more 
general nonholonomic systems, as well as applying them to quantize nonholonomic systems. 

2 Nonholonomic systems 

Nonholonomic mechanics takes place on a configuration space Q with a nonintegrable distribu- 
tion D that describes the (linear supposed) kinematic constraints of interest. These constraints 
are often given in terms of independent one-forms, whose vanishing in turn describes the dis- 
tribution. Moreover, one typically assumes that one can find a fibre bundle and an Ehresmann 
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connection A on that bundle such that D is given by the horizontal subbundle associated with 
A. Such an approach is taken, for example, in some recent books on nonholonomic systems 

PE]. 

Let Q be coordinatized by coordinates (r^, chosen in such a way that the projection of the 
above mentioned bundle structure is locally simply (r, s) ^ r. Moreover, let {w"} be a set of 
independent one-forms whose vanishing describes the constraints on the system. Locally, we can 
write them as 

a;"(r,s) = ds" + A?(r, s)dr^ 
The distribution D is then given by 

V = span{(9^/ — AJOgci}. 

One then derives the equations of motion using the Lagrange-d'Alembert principle, which takes 
into account the need for reaction forces that enforce the constraints throughout the motion of 
the system (see e.g. [3]). If L(r", s", r", s") is the Lagrangian of the system, these equations are 

d (dL\ dL _ d ( dL\ dL _ 

together with the constraints = —Afr^. One can easily eliminate the Lagrange multipliers A 
and rewrite the above equations in terms of the constrained Lagrangian 

Lc{r', s°,r^) = L{r^, -Afr^). 

The equations of motion, now in terms of Lc, become 

s° = -Afr^, 

d_(dLc\ ^^_j^a^_-j^a dL > (1) 



dt\dr^J dr^ ' ^ ds"" ds°^' 
where Bfj = d^jAf - d^iA'} + A'^jd.pA'} - A^jd.pAJ. 

To illustrate this formulation, consider perhaps the simplest example: a nonholonomically con- 
strained free particle with unit mass moving in M'^ (more details can be found in [3], In 
this example one has a free particle with Lagrangian and constraint given by 

L = i (x2 + y2 + i2) , i + xy = 0. (2) 

We can form the constrained Lagrangian by substituting the constraint into L, and proceed 
to compute the constrained equations, which take the form 

x = 0, y = -— — z = -xy. (3) 
1 + x"^ 

Another example of interest is the knife edge on a plane. It corresponds physically to a blade 
with mass m moving in the xy plane at an angle (p to the x-axis (see |24]). The Lagrangian and 
constraints for the system are: 

L =^m{x'^ +y'^) + ^J^'^, xsm{(l)) -ycos{(f>) =0, (4) 

from which we obtain the constrained equations: 

^ = 0, X = — tan{(j))(px, y = tan((/))x. 



3 



3 Second-order dynamics associated to a class of nonholonomic 
systems 



Recall from the introduction that we wish to investigate how we can associate a free Hamiltonian 
to a nonholonomic system. One way to do that is to rephrase the question in the Lagrangian 
formalism and to first investigate whether or not there exists a regular Lagrangian. Then, by 
means of the Legendre transformation, we can easily generate the sought after Hamiltonian. 
Rather than abstractly describing the various ways of associating a second-order system to a 
given nonholonomic system though, we will instead illustrate the method by means of one of 
the most interesting examples of a nonholonomic system. 

3.1 Associated Second-Order Systems for the vertically rolling disk 

The vertical rolling disk is a homogeneous disk rolling without slipping on a horizontal plane, 
with configuration space Q = M.'^ x x and parameterized by the coordinates {x,y,6,ip), 
where {x, y) is the position of the center of mass of the disk, 9 is the angle that a point fixed 
on the disk makes with respect to the vertical, and if is measured from the positive x-axis. The 
system has the Lagrangian and constraints given by 

X = Rcos{if)9, 

y = Rsm{^)9, (5) 

where m is the mass of the disk, R is its radius, and /, J are the moments of inertia about 
the axis perpendicular to the plane of the disk, and about the axis in the plane of the disk, 
respectively. The constrained equations of motion are simply: 

61' = 0, (^ = 0, x = Rcos{ip)9, y = Rsm{ip)9. (6) 

The solutions of the first two equations are of course 

9{t) = Ugt + 9o, (p{t) = u^t + (pQ, 

and in the case where u^p ^ 0, we get that the x- and y-solution is of the form 

x{t) = I — ) i?sin((/?(t)) + xo, 

y{t) = -(^^Rcos{ip{t))+yo, (7) 

from which we can conclude that the disk follows a circular path. If u^p = 0, we simply get the 
linear solutions 

x{t) = Rcos{Lpo)ugt + xq, y{t) = Rsm{(po)ugt + yo. (8) 

The situation in ([8]) corresponds to the case when remains constant, i.e. when the disk is 
rolling along a straight line. For much of what we will discuss in the next sections, we will 
exclude these type of solutions from our framework for reasons we discuss later. 
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Having introduced the vertical disk, let us take a closer look at the nonholonomic equations 
of motion ([6]). As a system of ordinary differential equations, these equations form a mixed 
set of coupled first- and second-order equations. It is well-known that these equations are 
never variational on their own O [7], in the sense that we can never find a regular Lagrangian 
whose (unconstrained) Euler-Lagrange equations are equivalent to the nonholonomic equations 
of motion ([1]) (although it may still be possible to find a singular Lagrangian). There are, 
however, infinitely many systems of second-order equations (only), whose solution set contains 
the solutions of the nonholonomic equations ([T|). We shall call these second-order systems 
associated second-order systems, and in the next section will wish to find out whether or not we 
can find a regular Lagrangian for one of those associated second-order systems. If so, we can use 
the Legendre transformation to get a full Hamiltonian system on the associated phase space. 
On the other hand, the Legendre transformation will also map the constraint distribution onto 
a constraint submanifold in phase space. The nonholonomic solutions, considered as particular 
solutions of the Hamiltonian system, will then all lie on that submanifold. 

There are infinitely many ways to arrive at an associated second-order system for a given non- 
holonomic system. We shall illustrating three choices below using the vertical rolling disk as an 
example. 

Consider, for example, taking the time derivative of the constraint equations, so that a solution of 
the nonholonomic system ([6]) also satisfies the following complete set of second-order differential 
equations in all variables {6,ip,x,y): 

9 = 0, (p = 0, X = —Rsm{{p)9ip, ij = Rcos{ip)6ip. (9) 

We shall call this associated second-order system the first associated second-order system. Ex- 
cluding for a moment the case where = 0, the solutions of equations ([9]) can be written 
as 

9{t) = uet + 9o 

ip{t) = U^t + (PQ 

^(*) = {^^^ Rsin{Lp{t)) + u^t + xq, 

y{t) = 



— Rcos{ip{tj) + Uyt + yo- 



By restricting the above solution set to those that also satisfy the constraints x = cos{ip)9 and 
y = sm{ip)9 (i.e. to those solutions above with Ux = Uy = 0), we get back the solutions ([7]) of 
the non-holonomic equations ([6]). A similar reasoning holds for the solutions of the form ([8]). 
The question we then wish to answer in the next section is whether the second-order equations 
([9]) are equivalent to the Euler-Lagrange equations of some regular Lagrangian or not. 

Now, taking note of the special structure of equations ([9|), we may use the constraints ([6]) to 
eliminate the 9 dependency. This yields another plausible choice for an associated system: 

" sin((/7) . . .. cos(ip) . . 

^ = 0, y5 = 0, x = f{x^, y = —^yip. (10) 

We shall refer to this choice later as the second associated second-order system. 

Lastly, we may simply note that, given that on the constraint manifold the relation sm{{p)x — 
cos{ip)y = is satisfied, we can easily add a multiple of this relation to some of the equations 
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above. One way of doing so leads to the system 

J(p = —mR{sm{ip)x — cos{ip)y)9, 

(I + mR^)9 = ■mR{sm{ip)x — cos{(p)y)(f, 

{I + mFf')x = —R{I + mR^) sm{(p)9(p + mR^ cos{(p){sm{ip)x — cos{(p)y)(p, 

{I + mR^)y = R{I + mR^) cos{(f)9ip + mR^ sm{(p){sm{ip)x — cos{Lp)y)(p. (11) 

For later discussion we shall refer to it as the third associated second-order system. We mention 
this particular second-order system here because it has been shown in \T6\ (using techniques 
that are different than those we will apply in this paper) that this complicated looking system 
is indeed variational! The Euler-Lagrange equations for the regular Lagrangian 

L = -^m(x^ + y2) + ^19'^ + ^ J^j^ + mR9{cos{(p)x + sm{ip)y), (12) 

are indeed equivalent to equations (|11|) . and, when restricted to the constraint distribution, its 
solutions are exactly those of the nonholonomic equations ([6]). We shall have more to say about 
this system in section 4.4 below. 



3.2 Associated Second-Order Systems in General 

We will, of course, not only be interested in the vertically rolling disk. It should be clear by now 
that there is no systematic way to catalogue the second-order systems that are associated to a 
nonholonomic system. If no regular Lagrangian exists for one associated system, it may still 
exist for one of the infinitely many other associated systems. For many nonholonomic systems, 
the search for a Lagrangian may therefore remain inconclusive. On the other hand, also the 
solution of the inverse problem of any given associated second-order system is too hard and too 
technical to tackle in the full generality of the set-up of the section 2. Instead, we aim here to 
concisely formulate our results for a well-chosen class of nonholonomic systems which include 
the aforementioned examples and for only a few choices of associated second-order systems. 

To be more precise, let us assume from now on that the configuration space Q is locally just 
the Euclidean space M" and that the base space of the fibre bundle is two dimensional, writing 
(ri,r2;Sa) for the coordinates. We will consider the class of nonholonomic systems where the 
Lagrangian is given by 

L=^ihrl + l2rl + ^Iasl), (13) 

a 

(with all Id positive constants) and where the constraints take the following special form 

Sa = -Aa{ri)r2. (14) 

Although this may seem to be a very thorough simplification, this interesting class of systems 
does include, for example, all the classical examples described above. We also remark that all 
of the above systems fall in the category of so-called Chaplygin systems (see f^). The case of 
2-dimensional distributions was also studied by Cartan, be it for other purposes (see e.g. [B] and 
the references therein). 
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In what follows, we will assume that none of the are constant (in that case the constraints 
are, of course, holonomic). The nonholonomic equations of motion ([T|) are now 

ri = 0, r2 = -N^{^If3Af3A'p)hr2, s„ = -^„r2, (15) 
where is shorthand for the function 

This function is directly related to the invariant measure of the system. Indeed, we have shown 
in [16] that for a two-degree of freedom system such as (|15p. we may compute the density 
of the invariant measure (if it exists) by integrating two first-order partial differential equations 
derived from the condition that the volume form be preserved along the nonholonomic flow. In 
the present case, these two equations read: 

and obviously the expression for N in (|16p is its solution up to an irrelevant multiplicative 
constant. In case of the free nonholonomic particle and the knife edge the invariant measure 
density is ~ l/\/T+~z? and ~ (1 -|- tan^((/))) = cos{(j)), respectively. In case of the 
vertically rolling disk it is a constant. We shall see later that systems with a constant invariant 
measure (or equivalently, with constant laA'^) always play a somehow special role. 

We are now in a position to generalize the associated second-order systems presented in section 
2.1 to the more general class of nonholonomic systems above. In the set-up above, the first 
associated second-order system is, for the more general systems ()15p . the system 

fi=0, r2 = -N^{^IpA,3A'p)hr2, = -{A'^rir2 + Aar2), 

or equivalently, in normal form, 

ri = 0, ra = -AT^ ( ^ I^A^4)rir2, 

A'^-N^A^Y^IpApA'^)) rir2. (18) 

P 

For convenience, we will often simply write 

fi =0, r2= T2{ri)fif2, Sa = ^airi)rir2, 
for these types of second-order systems. 

The second associated second-order system we encountered for the vertically rolling disk also 
translates to the more general setting. We get 

n = 0, fs = -N'^{^I,3ApA'p)rir2, 

ya-N'A4^I,A,A',))n{^^y (19) 
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where in the right-hand side of the last equation, there is no sum over a. A convenient byproduct 
of this way of associating a second-order system to (fT5|) is that now all equations decouple except 
for the coupling with the ri-equation. To highlight this, we will write this system as 

r\ = 0, Qa = Ha(ri)gan 

(no sum over a) where, from now on, (qa) = {r2,Sa) and (qi) = (ri,^^). 

We postpone the discussion about the third associated second-order system of our class until 
section 4.4. 



4 Lagrangians for associated second-order systems 
4.1 The inverse problem of Lagrangian mechanics 

Let Q be a manifold with local coordinates (g*) and assume we are given a system of second- 
order ordinary differential equations g* = f^{q,q) on Q. The search for a regular Lagrangian 
is known in the literature as 'the inverse problem of the calculus of variations,' and has a long 
history (for a recent survey on this history, see e.g. [21] and the long list of references therein). 
In order for a regular Lagrangian L{q,q) to exist we must be able to find functions gij{q,q), 
so-called multipliers, such that 

It can be shown [21 \T2\ [27] that the multipliers must satisfy 

det{gij) / 0, 

where V*- = — ^S^j/* and 

= r (^d^,f) - 2d^,f - \d^,fd^,f\ 

The symbol V stands for the vector field q^dqi + f^d^i on TQ that can naturally be associated 
to the system = P{q,q). Conversely, if one can find functions gij satisfying these conditions 
then the equations = /* are derivable from a regular Lagrangian. Moreover, if a regular 

Lagrangian L can be found, then its Hessian . is a multiplier. 

oq^oq^ 

The above conditions are generally referred to as the Helmholtz conditions. We will fix from the 
start gtj = gji for j <i, and we will simply write gijk for d^kgij, and also assume the notation 
to be symmetric over all its indices. 

The Helmholtz conditions are a mixed set of coupled algebraic and PDE conditions in (gij). 
We will refer to the penultimate condition as the 'V- condition,' and to the last one as the 
'^-condition.' The algebraic ^-conditions are of course the most interesting to start from. In 



gji — 9ij 1 



dgij ^ dg_ 
dq^ dqi 



'ik 
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fact, we can easily derive more algebraic conditions (see e.g. [llj). For example, by taking a F- 
derivative of the ^-condition, and by replacing T{gij) everywhere by means of the V-condition, 
we arrive at a new algebraic condition of the form 

gik{V^)j = gjkiv^)^ 

where (V$)*- = r($*) - V^<i>^ - V^$j„. As in ^Jjj, we will call this new condition the 
(V^>)-condition. It will, of course, only give new information as long as it is independent from 
the ^-condition (this will not be the case, for example, if the commutator of matrices [<I>, V<&] 
vanishes). One can repeat the above process on the (V<l>)-condition, and so on to obtain possibly 
independent (V . . . V<I>)-conditions. 

A second route to additional algebraic conditions arises from the derivatives of the ^-equation 
in (^-directions. One can sum up those derived relations in such a way that the terms in gijj- 
disappear on account of the symmetry in all their indices. The new algebraic relation in gij is 
then of the form 

dtjKi + aijRlk + dkjRji = 0, 

where Rj,i = dgj{^^) — dqi{^j). For future use, we will call this the i?-condition. 

As before, this process can be continued to obtain more algebraic conditions. Also, any mixture 
of the above mentioned two processes leads to possibly new and independent algebraic conditions. 
Once we have used up all the information that we can obtain from this infinite series of algebraic 
conditions, we can start looking at the partial differential equations in the V-conditions. 

We are now in a position to investigate whether a Lagrangian exists for the two choices of 
associated systems ([T8]) and (fT9|) . 

4.2 Lagrangians for the first associated second-order system 

The first second-order system of interest is of the form 

h = 0, r2 = F2(ri)rir2, Sa = Fa(ri)rir2- (20) 
The only non-zero components of (<!>*•) are 

^? = (^r2-F'2)nr2, ^l = -{lrl-r',)fl 

CD? = (iF„F2 - K)hr2, '^^ = -(^r,r2 - T'Jrl 
For V<I> and VV<I> we get 

(v$)? = (r2r'2 - r'2')r?r2, (v<i>)i = -(F2F'2 - r'2')r?, 
(v$)? = (r„r2 - r:;)r?r2, (v#)f = -(r,r2 - r'yi 

and 

(VVd))? = ((r'2)2 + F2F^'-F^>?r2, 

{vv^)l = _((r'2)2 + r2r'2'-rr)rt 

(vv<i>)f = (r'„r'2 + ^r,r(,'-ir::r2-r>?r2, 

(VV$)f = -(r:,r'2 + ^F„F'2'-iF'^F2-F-)rt, 
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and so on. 



We can already draw some immediate consequences just by looking at the above explicit expres- 
sions. Let's make things a bit more accessible by considering the case where the dimension is 4. 
Then, the ^-equations of the system (j20p and their derivatives are all of the form 

5l2*i + fflS^i + 914^2 = 522*? + 523*? + 524*1, 
523*? +533*?+ 534*? = 0, 

523*1 + 533*i + 534*1 = 0, (21) 
524*? + 534*? + 544*? = 0, 
524*1 + 534*i + 544*2 = 0, 

where, within the same equation, ^ stands for either V$, VV<I>, VVV$, ... We will refer 
to the equations of the first line in ([2T]) as 'equations of the first type,' and to equations of the 
next four lines as 'equations of the second type.' The first 3 equations of the first type, namely 
those for V$ and VV<I> are explicitly: 

512^1 + 5l3*i + 514^2 = 522^? + 523^? + 524^?, 

5i2(V*)l + 5i3(V^>)i + 5l4(V$)t = 522(V«>)? + 523(V«')? + 524(V$)?, (22) 
ffl2(VV$)l + 5l3(VV$)i + 514(VV^>)^ = (722(VV$)? + 523 (VV$)? + 524(VV«>)?. 

For the systems at hand, the particular expression of $ and its derivatives are such that 

$1(V$)? -$?(V$)1 = 0, 
(V«>)1(VV«>)? - (V^>)?(VV^>)1 = 0, 

and so on. By taking the appropriate linear combination of the first and the second, and of 
the second and the third equation in (j22p . we can therefore obtain two equations in which the 
unknowns gi2 and g'22 are eliminated. Moreover, under certain regularity conditions, these two 
equations can be solved for gi^ and 514 in terms of 523 and 524 (we will deal with exceptions 
later on). So, if we can show that (723 and 524 both vanish, then so will also 513 and 514. Then, 
in that case 512*2 = 522*?, but no further relation between 512 and 522 can be derived from 
this type of algebraic conditions. 

The infinite series of equations given by those of the second type in (j2ip are all equations in 
the 5 unknowns 523, 533 , 534 , 524 and 544. Not all of these equations are linearly independent, 
however. In fact, given that the system ()20p exhibits the property 

(where ^ is one of V<I>, VV<I>, ...), one can easily deduce that the last four lines of equations 
in (j2ip actually reduce to only two kinds of equations. If we assume that we can find among 
this infinite set 5 linearly independent equations, there will only be the zero solution 

523 = 533 = 534 = 524 = 544 = 0, 

and from the previous paragraph we know that then also 514 = gi3 = 0. To conclude, under the 
above mentioned assumptions, the matrix of multipliers 





511 


512 









512 


522 






















I 











/ 



10 



is singular and we conclude that there is no regular Lagrangian for the system. The above 
reasoning can, of course, be generalized to lower and higher dimensions. 

We will refer to the above as 'the general case'. The assumptions made above are, however, 
not always satisfied, and they need to be checked for every particular example. Let us consider 
first the example of the (three-dimensional) nonholonomic particle, where r2 = —x/il + x^) and 
Fa = — 1/(1 + 2;^). The equations for ^ = <1>, V<1> of the second type give the following two linear 
independent equations 

(x^ - 2)523 + 3x533 = 0, (x^ - 5x)523 + (5x2 _ -^^^^^ ^ Q 

We can easily conclude that 523 = 533 = 0. With that, the first two equations of the first type 
are 

(x^ - 2)xgi2 + 3xx5ri3 + (x^ - 2)1/522 = 0, 

(x^ - 5x)x5i2 + (5x2 _ i^±g^^ (2;3 _ 5x)y522- 

From this 513 = and X512 = — 1/922, and there is therefore no regular Lagrangian. 

With a similar reasoning (but with different coefficients) we reach the same conclusion for the 
example of the knife edge on a plane. 

The vertically rolling disk is a special case, however, and so is any system (llSp with the property 
that Y2a -^cf^a is * constant. This last relation is in fact equivalent with the geometric assumption 
that the density of the invariant measure is constant. In that case, we get r2 = 0. Not only 
does r2 vanish, but so do all ^'f and ^'2 ^ — V<I>, .... We also have r3 = —Rsin{ip) and 
r4 = Rcos{ip). Moreover, looking again first at expressions (j22p . one can easily show that for 
the vertically rolling disk these three equations, and in fact any of the equations that follow in 
that series, are all linearly depending on the following two equations 

cos(v9)05i3 + sin((^)05i4 + cos{ip)0g23 + sin((^)^524 = 0, 
sm{ip)Lpgi3 - cos{ip)ipgii + sin(v9)^523 - cos((/?)^524 = 0. 

Although these equations are already in a form where 512 and 522 do not show up, it is quite 
inconvenient that there is no way to relate these two unknowns to any of the other unknowns. 
However, as in the general case, we can deduce from this an expression for 513 and 514 as a 
function of 523 and 524- We get 

9 
gi3 = — r523, 914 = — r524- (23) 

The infinite series of equations of the second type (i.e. the last four lines in (j21|) ) are all linearly 
dependent to either one of the following four equations 

COs(<y9)533 + sin(99)534 = 0, COs{ip)g3A + Sin(v9)544 = 

sin(v')533 - cos(93)534 = 0, sm{(p)g34 - cos{(p)g44 = 0, 

from which 533 = 534 = 544 = follows immediately. In comparison to the general case, however, 
we can no longer conclude from the above that also 523 and 524 vanish, and therefore, we can 
also not conclude from (j23p that 513 and 514 vanish. This concludes, in fact, the information 
we can extract from the <I>-condition, and the algebraic conditions that follow from taking its 
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derivatives w.r.t. V. Also, any attempt to create new algebraic conditions by means of the tensor 
R is fruitless, since an easy calculation shows that, when the above conclusions are taken already 
into account, all equations that can be derived from R are automatically satisfied. However, 
we have enough information to conclude that there does not exist a regular Lagrangian for the 
vertically rolling disk and its first associated second-order system. Indeed, the determinant of 
the multiplier matrix 

5^12 
922 
923 
92A 



( 



i9ij) 



911 

912 
Afi'23 



^923 
923 






A524 

524 






(with A = —9/ip) clearly vanishes and this is a violation of one of the first Helmholtz conditions. 

Thus, to summarize the above results, for the nonholonomic free particle the knife edge on 
the plane (|4|) and the vertically rolling disk ([5]), we conclude that there does not exist a regular 
Lagrangian for their first associated second-order system (jlSp . 



4.3 Lagrangians for the second associated second-order system 

In this section, we will investigate the inverse problem for the second associated system, 

n = 0, Qa = ^a{ri)qari- (24) 

In the (^a-equations, there is no sum over a, which is an index that runs from 2 to the dimension of 
the configuration space, and with respect to the formulation of the inverse problem in section 3, 
we have /i = and fa = EaQaTi- Moreover, one can easily compute that the only non- vanishing 
components of ^ are now 

The ^-conditions turn out to be quite simple: if <^J?[ 7^ 0, then 

qa9aa = -ri9la, (25) 

and if / ^1 « / then 

9ab = 0. (26) 

These restrictions on $ lead to the assumptions that first Hq 7^ and 7^ 2/(C — ri), where 
C is any constant, second that 7^ Ef, and, formally, — H;, 7^ E}j/{C — j Ehdri), where 
Ebi^i) = exp(J 2Ei,dri). Suppose for now that we are dealing with nonholonomic systems ([T9]) 
where this is the case. Then one can easily show that all the other V . . . V<I>-conditions do 
not contribute any new information, as well as that the -R-condition is automatically satisfied. 
Thus we should therefore turn our attention to the V-condition, which is a PDE. To simplify 
the subsequent analysis though, we note that although the multipliers gij can in general be 
functions of all variables {ri,qa, riiila)-, in view of the symmetry of the system we shall assume 
them to be, without loss of generality, functions of (ri,ri,ga) only. 

Now, by differentiating the algebraic conditions by ri, ri and qa, we get the additional conditions 

qa9'aa = -^l9la 

9aa + Qa9aaa = —ri9laa, Qa9laa = — 9la — ^igua 
9aab = = 9lab, H a ^ b. 
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Finally, the V-Helmholtz conditions are, with the above already incorporated, 

g'li + X] '^bi.giibQb - gbb^) = o, 
b ^1 

9aa + '^a{gaaa<ia + Qaa) = 0. 

In what follows we will implicitly assume everywhere that ri 7^ 0. As a consequence, the 
multipliers (gij) (and the Lagrangians we may derive from it) will only be defined for ri 7^ 

It is quite impossible to find the most general solution for {gij) though. We will show that there 
is an interesting class of solutions if we make the anszatz that gf^f,}, = for all b. With that and 
with the above gaab = in mind, we conclude that all such gbi, will depend only on possibly ri 
and fi. Moreover, from the last V-conditions we can determine their dependency on the variable 
ri. Since now 

9'bb + dbb'^b = 0, 

it follows that gbb{fi-,'i'i) = -^^(n) exp(— ^5(ri)), where is such that .^^ = H5 and where Fb{ri) 
is still to be determined from the remaining conditions. From one of the above conditions we 
get gibb = -gbb/h (since gm = 0), so 

dfi f\ ' 

from which = C^/ri, with Cf, a constant, and thus gi,i, = exp(— ^ft)/ri. Therefore, from the 
algebraic conditions, = -igbb/h)qb = -Cbew{-Cb)Qb/ri, and thus = 2C6^6exp(-^fe)/rf . 
With this, the first V-condition becomes 

b ''1 

and thus 

gii = ^Cbe-x.^{~^,b)^ + C{ri,qb)- 
b ""i 

Given that gn^ = 2C6gbexp(— ^b)/rf, we can now determine the gfe-dependence of C. We simply 
get ^ 

gii = Y,Cb^M-^b)% + Fi{h). 

b '^i 

Notice that gm does not show up explicitly in the conditions or in the derived conditions. 
Therefore, there will always be some freedom in the ^n-part of the Hessian, represented here 
by the undetermined function Fi (ri ) . 



Up to a total time derivative, the most general Lagrangian whose Hessian gij = ^ ..^ ,. is the 
above multiplier, is: 



L = Pih) + IY1 CbeM-^b)f, (27) 
^ b ^1 



where cfp/drf = Fi. One can easily check that the Lagrangian is regular, as long as cfp/drf 
is not zero, and as long as none of the Cb are zero. Remark, finally, that the Lagrangian is 
only defined on the whole tangent space if = (and p is at least everywhere). We can 



13 



therefore only conclude that there is a regular Lagrangian (with the ansatz (^tfcb = 0) on that part 
of the tangent manifold where ri 7^ 0. As a consequence, the solution set of the Euler-Lagrange 
equations of the Lagrangian ()27p will not include those solutions of the second-order system (I20p 
where ri = 0. In case of the vertically rolling disk, for example, these solutions are exactly the 
special ones given by ([7]), and that is the reason why we will exclude them from our formalism. 

Recall that at the beginning of this section, we have made the assumptions that ^ and 
/ ^b- Suppose now that one of these assumptions is not valid, say H2 = and therefore 
$^ = 0. Then, among the algebraic Helmholtz conditions there will no longer be a relation in 
(j25p that links §22 to §12- In fact, since the gij now need to satisfy only a smaller number of 
algebraic conditions, the set of possible Lagrangians may be larger. We can, of course, still take 
the relation 

92522 = -hgu (28) 

as an extra ansatz (rather than as a condition) and see whether there exists Lagrangians with 
that property. By following the same reasoning as before, we easily conclude that the function 
(j27p is also a Lagrangian for systems with <I>| = 0. In fact, it will be a Lagrangian if any of the 
assumptions is not valid. 

Apart from ([28]) . we are, of course, free to take any other ansatz on gi2 and §22- If we simply 
set 

512 = 0, 



it can easily be verified that also 



1 . . 



L = p{h)+a{r2) + -\^}^a^exp{-^^)^j (29) 

is a Lagrangian for a system (p^ with H2 = (where, as usual, (qa) = (^2, Sq,)). It is regular as 
long as both d?p/drf and d?a/dr2 do not vanish. 

Proposition 1. The function 

L = ,in) + ^[cS^Ec,^^. (30) 

with d?p/drf ^ and all Ca ^ is in any case a regular Lagrangian for the second associated 
systems il9\). If the invariant measure density N is a constant, then also 

L = piri) + air2) + ^Y.^,^^, (31) 

where d'^p/drf / 0, d^a/dr\ / and all Ca ^ is a regular Lagrangian for the second associated 
systems ( f7g|) . 



Proof. For the second associated systems, the second-order equations (|24p are of the form (|19p . 
One easily verifies that in that case 

^2 = lniV and C„ = ln(iVylc.) (32) 

are such that = H^. The first Lagrangian in the theorem is then equal to the one in (j27p . For 
a system with constant invariant measure N , we get that H2 = 0. Therefore, also the function 
([29l) is a valid Lagrangian. □ 
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Let us end this section with a hst of the Lagrangians for the nonholonomic free particle, the 
knife edge on a horizontal plane and the vertically rolling disk. The respective Lagrangians pop 
for the first two examples are: 

L = p{x) + \^Jl + x^ (C2^ + C3—) , (33) 

and 

L = pi^) + l^m{l + tan(0)2) (c^t- + ^3-^) , 
^ V <P tan((p)(p/ 



p(<A) + ^C^2^/^^^ + ICsV^-y^. (34) 
^ cos ^ smi 



The vertically rolling disk in one of those systems with constant invariant measure. The first 
Lagrangian (f30]) is: 



,-,2 



L = ,(^) + :^Z+™«! C.«; + C3^ + C ' ) (35) 

2 \ cos(99j(^ sm[ip)(p I 



and the second Lagrangian ([3T]) is: 



L = ;.(vi) + a(^) - ^^I^^ f «3^^ + a4-4V ) • (36) 

2 \ cos((/jjv9 sm(99j99/ 



4.4 Lagrangians for the third associated second-order system 

In section 2.1 we have described a third associated second-order system (|lip for the example 
of the vertically rolling disk. That system comes actually from a comparison of the variational 
nonholonomic and the Lagrange-d'Alembert nonholonomic equations of motion we conducted in 
|16j . There we investigated the conditions under which the variational nonholonomic Lagrangian 
Ly would reproduce the nonholonomic equations of motion when restricted to the nonholonomic 
constraint manifold. Thus, instead of associating second-order systems to nonholonomic equa- 
tions and applying the techniques of the inverse problem to derive the Lagrangian (and the 
Hamiltonian), in ^6] we started from a specific Lagrangian (the variational nonholonomic La- 
grangian Ly) and investigated the conditions under which its variational equations match the 
nonholonomic equations. Other relevant work on this matter can be found in e.g. [8]. 

In case of our class of nonholonomic systems with Lagrangian ()13p and constraints (jl4p the 
variational nonholonomic Lagrangian is simply 

Lv = L -^^-^{80 + Aar2) 

= -{Iirl + hrl -^^lasi) -^^Ao^IaSah. 

a a 

A short calculation shows that its Euler-Lagrange equations in normal form are given by 

P p p 

Sa = -(< - A^2^4 j;//3^;34))rir2 -^a(^//34i/3)ri. (37) 

p p 
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In general, these systems are not associated to our class of nonholonomic systems. That is, the 
restriction of their solutions to the constraint manifold Sa = —■A.a'f'2 are not necessarily solutions 
of the nonholonomic equations (llSp . However, in case that the invariant measure density is a 
constant, we have that Ylp-^p-^p-^'p — As a consequence, all the terms in the equations ([37|) 
that contain IpA'^^sp vanish when we restrict those equations to the constraint manifold and 
the equations in Sa integrate to the equations of constraint (jl4p . The restriction of the equations 
(j37p is therefore equivalent with the nonholonomic equations (115^ . We conclude the following. 

Proposition 2. // N is constant, the equations JgTp form an associated second-order system 
and, by construction, they are equivalent to the Euler- Lagrange equations of the variational 
nonholonomic Lagrangian Ly- 

We refer to [16] for more details and some more general statements on this way of finding 
a Lagrangian for a nonholonomic system and we end the discussion on the third associated 
systems here. 



5 Hamiltonian formulation and the constraints in phase space 

In the situations where we have found a regular Lagrangian, the Legendre transformation leads 
to an associated Hamiltonian system. Since the base solutions of the Euler-Lagrange equations 
of a regular Lagrangian are also base solutions of Hamilton's equations of the corresponding 
Hamiltonian, the Legendre transformation FL will map those solutions of the Euler-Lagrange 
equations that lie in the constraint distribution T) to solutions of the Hamilton equations that 
belong to the constraint manifold C = FL{T>) in phase space. Recall however that the La- 
grangians for the second associated second-order systems (and their Legendre transformation) 
were not defined on ri = 0, and so will also the corresponding Hamiltonians. 

Let us put for convenience p{ri) = \liri and cr(r2) = \l2'f'2 ™- the Lagrangians of the previous 
section. 

Proposition 3. Given the second associated second-order system il9\) . the regular Lagrangian 
llcl(J\) ( away from ri = Q) and constraints j j^p on TQ are mapped by the Legendre transform to 
the Hamiltonian and constraints in T*Q given by: 



(38) 



In case N is constant, the second Lagrangian [31\) and constraints p4\ ) are transformed into 

H = ^^pI + ^^{pi + \^ ^pI^ j ' I^^^iVa + aaP2 = 0, (39) 

where riin,pi,Pa) = {Pl + ^^J2a'^aPa/'^a)/h- 

Proof. The Legendre transformation gives for the Lagrangian ()27p 

Pl = /in - ^ V'Cfcexp(-^fc)%, ph = Cbexp{-^b)^, (40) 
2 ^-^ rf ri 
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from which one can easily verify that the corresponding Hamiltonian is 

H = ^\pi + ;:} eMd,)^] ■ (41) 



In the case of the second associated second-order systems in the form (|19p . the take the form 
(j32p . and we obtain the Hamiltonian in expression (j38p . From (j40p we can then compute the 
constraint manifold C in phase space. Since now 

P2 = (-^2^rn- and Pq, = Cq 



■iVri "iVri' 

the constraints ()14p can be rewritten as 

where ri = 7|-(pi + \N{j)2r\IC2 + ^/3P^/C'/3))- Assuming as always that r\ / 0, we get that 
the constraint manifold in phase space is given by C^Va = —CaP2 for all a. 



An analogous calculation with the Lagrangian (j3ip gives the Hamiltonian and the constraints 
in ()39p . in the case where is constant. □ 



We can recover the Hamiltonians of [T7j from Proposition 2. As perhaps the simplest example, 
note that with {ri,r2,Sa) = {x,y,z), by taking C2 and C3 both to be 1, and A{ri) = x, 
we recover the Hamiltonian and the constraint that appears in [17] for the nonholonomic free 
particle. 

Consider now the knife edge on the plane. Taking C2 = C3 = l/^/m and A{(j)) = — tan((/)) gives: 

H = ^(^Pc^ + l{cos{^)pl-sm{(P)Py?j , (42) 
while the constraint manifold becomes 

Px+Py = 0- (43) 
For the rolling disk we get for the first Hamiltonian (j38p 

p1 cos{lp)pI sin(v?)p2 





2/rTT^ \ C2 C3 C4 



and C2PX = —CsPe and C2Py = —C/^pq for the constraints. These are not the Hamiltonian 
and the constraints that appear in ^17j though. It turns out that the Hamiltonian and the 
constraints in [ITj are in fact those that are associated to the second Hamiltonian (j39p . It is, 
with, for example, 03 = 04 = — J/\// + mE? of the form 

1 1/1 1 \ ^ 

H = —pj + 2 ( + ^°^('^) + 2^y ^^'^('^) ) 

and the constraints are 

V'Px=Pe, •pPy=Pe 
where <p = Pip + \ cos{ip)p'^ + ^ sm{ip)py or, equivalently, 

Px-Py = 0, ippx -pe = 0, 

as the constraints appears in pLZj • 
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6 Pontryagin's Maximum Principle 

Consider the optimal control problem of finding the controls u that minimize a given cost function 
G{x, u) under the constraint of a first order controlled system x = f{x, u). One of the hallmarks 
of continuous optimal control problems is that, under certain regularity assumptions, the optimal 
Hamiltonian can be found by applying the Pontryagin maximum principle. Moreover, in most 
cases of physical interest, the problem can be rephrased so as to be solved by using Lagrange 
multipliers p. Form the Hamiltonian H^{x,p,u) = {p,f{x,u)) — poG{x,u) and calculate, if 
possible, the function u*{x,p) that satisfies the optimality conditions 



Then, an extremal x{t) of the optimal control problem is also a base solution of Hamilton's 
equations for the optimal Hamiltonian given by H*{x,p) = H^{x,p,u*{x,p)). The optimal 
controls u*{t) then follow from substituting the solutions {x{t),p(t)) of Hamilton's equations for 
H* into u*{x, q). 

Such a usage of the multiplier approach can also be applied with succes to the mechanics of phys- 
ical systems with holonomic constraints. However, in the case of nonholonomically constrained 
systems the Lagrange multiplier approach, also called the vakonomic approach by Arnold [l], 
generally leads to dynamics that do not reproduce the physical equations of motion (see [HI [22] 
and references therein). Thus, the rich interplay between Pontryagin's Maximum Principle, the 
vakonomic approach, and the physical equations of motion of a constrained system breaks down 
when the constraints are nonholonomic. However as we showed in a previous paper [16], for 
certain systems and initial data the vakonomic approach and Lagrange-D'Alembert principle 
yield equivalent equations of motion. 

We will show here for the second associated systems 

n = 0, qa = '^a{ri)qari, 

that we can also find the Hamiltonians of the previous section via a rather ad hoc application 
of Pontryagin's Maximum Principle. Hereto, let us put = as before and observe that the 
above second-order system can easily be solved for {ri{t),qa{t)). Indeed, obviously ri is constant 
along solutions, say ui. We will suppose as before that ui / 0. From the ga-equations it also 
follows that (ja/ exp(^a) is constant, and we will denote this constant by Ua- To conclude, 

hit) = Ul, qa{t) = Uaexp(^a(ri(t))). 

Keeping that in mind, we can consider the following associated controlled first-order system 



(no sum over a), where (ui,Ua) are now interpreted as controls. 

The next proposition relates the Hamiltonians of Proposition 2 to the optimal Hamiltonians for 
the optimal control problem of certain cost functions, subject to the constraints given by the 
controlled system (I44p . 

Proposition 4. The optimal Hamiltonian H* of the optimal control problem of minimizing the 
cost function 



{x,p, u*{x,p)) = 0. 



du 



h = Ul, qa = Ua exp(^a(ri)) 



(44) 
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subject to the dynamics ( (^^[ j is given by: 



(45) 



If H2 is zero, the optimal Hamiltonian for the optimal control problem of minimizing the cost 
function 

subject to the dynamics 1^44^ ^-^ given by: 



(46) 



In case the controlled system is associated to a nonholonomic system (that is, in case the 
take the form ^3^), the above Hamiltonians are respectively the Hamiltonians [38\) and [W\) of 
Proposition 2. 



Proof. The Hamiltonian is 

H^{ri,qa,Pl,Pa,Ui,Ua) = PlUi + y^PaUaexp(^a) - Gi- 



(47) 



The optimahty conditions dH^/dui = 0, dH^/dua = 0, together with the assumption that 
ui 7^ 0, yield the following optimal controls as functions of {q,p): 



hul = pi + ij^exp(ea)^ 



Pa_ 
Ca 



For the Hamiltonian H*{q,p) = H^{q,p,u*{q,p)), we get 
H*{q,p) = 



1 

z 



Pi - \hu{^ ul + ^exp(e„)< (^pa - ^Ca^^ 



which is exactly the Hamiltonian (j4ip . 

For the second cost function, with H2 = 0, we get for Pontryagin's Hamiltonian 

= PnUri +Pr2Ur2 + ^ Pq exp (^^ ) - G2. 
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1 

^ / ^ — ' n ^ 



The optimal controls as functions of {q-,p) are now 

a 

K _ Pa 

With this the Hamiltonian becomes 



which is exactly (|39|) after the substitution (j32|) . □ 



7 Related Research Directions and Conclusions 

In essence, the method we have introduced in the previous sections resulted in an unconstrained, 
variational system which when restricted to an appropriate submanifold reproduces the dynamics 
of the underlying nonholonomic system. Although we have restricted our attention to a certain 
explicit subclass of nonholonomic systems, many of the more geometric aspects of the introduced 
method seem to open the door to generalizing the results to larger classes of systems. For 
example, a lot of what has been discussed was in fact related to the somehow hidden symmetry 
of the system. That is to say: both the Lagrangian ()13p and the constraints (I14p of the systems 
at hand were explicitly independent of the coordinates r2 and Sa- This property facilitated the 
reasoning we have applied in our study of the corresponding inverse problems. One possible path 
to the extension of some of the results in this paper may be the consideration of systems with 
more general (possibly non-Abelian) symmetry groups. A recent study [10] of the reduction of 
the invariant inverse problem for invariant Lagrangians may be helpful in that respect. 

The methods of the inverse problem have lead us to the Lagrangians for the second associated 
second-order systems ()24p . For those systems the g"-equations were, apart from the coupling 
with the ri-equation, all decoupled from each other. It would be of interest to see, for more 
general systems, how such a form of partial decoupling influences the question of whether or not 
a regular Lagrangian exists. 

In the previous section we have found a new link between the fields of optimal control, where 
equations are derived from a Hamiltonian, and nonholonomic mechanics, where equations are 
derived from a Hamiltonian and constraint reaction forces. By combining elements of both 
derivations, for certain systems one can formulate the mechanics in a form analogous to the 
treatment of constraints arising from singular Lagrangians that leads to the Dirac theory of 
constraints [18] , which allows for the quantization of constrained systems wherein the constraints 
typically arise from a singular Lagrangian (see |20j and references therein). Central to the 
method is the modification of the Hamiltonian to incorporate so-called first and second class 
constraints. 

The method proposed in this paper in a sense provides an analogue to Dirac's theory and 
allows for the investigation into the quantization of certain nonholonomic systems by similarly 
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modifying the usual Hamiltonian. In attempting to quantize the class of systems we have 
considered, we can now instead use one of the Hamiltonians found in Proposition 3. We should 
note that there have already been some attempts to quantize nonholonomic systems [U [T3| \T7\ 
[T9 \ 120 1 [25]. f-iid that the results have been mixed, mainly due to the inherent difficulties arising 
in the quantization procedure, as well as the difficulties in dealing with the system's constraints. 
However, the present work enables one to treat the constraints more like an initial condition, 
since, for example, the constraint ()43p is really the relation ci + C2 = 0, where px = ci, and 
Py = C2 follows from H in (|42p being cyclic in x, y. Such a treatment of the constraints eliminates 
much of the difficulty arising in attempting to quantize some nonholonomic systems. 

As example, consider the knife edge on the plane, in view of (j42p . We can take the quantum 
Hamiltonian H to be of the form 



which is a Hermitian operator, and consider the quantum version of the constraint manifold (|43p 
Px Py = (Here " stands for the quantum operator form under canonical quantization) . There 
have in the literature been essentially two different ways to impose these quantum constraints: 
strongly and weakly. One may require that the quantum constraints hold strongly, by restricting 
the set of possible eigenstates of ()48p to those which satisfy the quantum operator form of 
constraint ()43p . On the other hand, one may only require that the eigenstates \^) of (|48p satisfy 
the quantum constraints in mean 



a weaker condition but arguably a more physically relevant viewpoint also advocated in |17tll9j. 
In [T7| the authors show that (using the example of the vertical rolling disk) the weaker version 
can be used to recover the classical nonholonomic motion in the semi-classical limit of the 
quantum dynamics. Details of our application of the methods in this paper to the quantization 
of nonholonomic systems will be presented in a future publication. 
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